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Abstract 

The spatial string tension across a crossover from the low temperature phase to the high temper- 
ature phase is computed in QCD with two flavors of non-perturbatively improved Wilson fermions 
at small lattice spacing a ~ 0.12 fm. We find that in the low temperature phase spatial string 
tension agrees well with zero temperature string tension. Furthermore, it does not show increasing 
for temperatures up to T = 1.36T pc , the highest temperature considered. Our results agree with 
some theoretical predictions. 

PACS numbers: ll.15.Ha, 12.38.Aw, 12.38.Gc 
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I. INTRODUCTION 



In view of running and future heavy-ion collisions experiments determination of the 
(pseudo) critical temperature T pc as well as studies of other features of the finite temper- 
ature transition in full QCD with Nf = 2 and Nf = 2 + 1 flavors of dynamical quarks are 
subjects of many recent lattice investigations. Three groups employed improved staggered 
fermions and an improved gauge field action |l|, 2 , [J . Simulations with staggered fermions, 
although computationally less demanding, are made with the use of the fourth root of the 
fermion determinant. This introduces theoretical uncertainty. Thus studies with the Wilson- 
type quark actions free of such problem are necessary. The CP-PACS 4], WHOT-QCD \^ 
and DIK 6| collaborations used improved Wilson fermions with an improved or unimproved 
gauge field action. In this paper we are using the lattice configurations generated by DIK 
collaboration. 

Recent results jl, 6] show importance of the finite lattice spacing effects in the studies of 
the finite temperature transition in QCD. The DIK collaboration simulations were done on 
lattices with large N t , i.e. the transition was studied at small lattice spacing, close to the 
continuum limit. For this reason we consider it important to study in more details various 
features of the finite temperature QCD using gauge field configurations accumulated by DIK 
collaboration. In this letter we present results for the spatial string tension a s . 

In the presence of dynamical quarks the flux tube formed between static quark- ant iquark 
pairs breaks at large distances. In contrast the so called spatial Wilson loop, i.e. the Wilson 
loop with both directions being spatial is found to show non-vanishing spatial string ten- 
sion even at T > T pc . This is called the spatial confinement. a s is one of the fundamental 
quantities to characterize the perturbative and non-perturbative properties of the hot QCD 
medium. It provides useful information about the effective 3D model, which can be obtained 
by integrating out heavy modes. It is believed that non-perturbative contributions to ob- 
servables, which at high temperature are dominated by contributions of magnetic modes, 
are flavor independent apart from indirect dependence through the flavor dependence of the 
running coupling. a s can be used to test this aspect of the dimensional reduction approach 
to QCD. 

Although quarks are expected to decouple from the spatial observables for T ^> T pc due 
to dimensional reduction and thus do not affect a s in the high temperature limit, it is not 
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TABLE I: The parameters of lattice configurations used in computations of a s . is the physical 

string tension dependent on m^ro, respective values are taken from [lij] 

obvious whether the same is true near T pc . Thus lattice computations are also useful to 
check theoretical predictions for this range of temperatures. 



II. DETAILS OF COMPUTATIONS 



o s was recently computed by RBC-Bielefeld 17 
Our lattice action is different from that used in [7, 



8j and WHOT-QCD [5) collaborations, 
and similar to action used in [511 . The 



important difference from jlj] is that we are working at two times smaller lattice spacing, i.e. 
much closer to the continuum limit. We perform computations on N% N t = 16 3 8 lattices 
while in j^] N t = 4 lattices were used. The configurations were generated with the usual 
Wilson gauge field action and the non-perturbatively 0(a) improved Wilson fermionic action 
S F : 

S F = Sf ] - ~k g c sw a 5 V ]ip(s)a tlu F flu (s)ilj(s) } (1) 



where Sp is the original Wilson fermionic action, g is the gauge coupling and F^ v (x) is the 
field strength tensor. The clover coefficient c sw is determined non-perturbatively [9(. We use 
configurations generated at two values of the lattice coupling f3 = Q/g 2 : j3 = 5.2 and (3 = 
5.25, for temperatures between 0.78 T pc and 1.36 T pc . To fix the physical scale and the quark 



mass results obtained with this action by QCDSF and UKQCD collaborations 



3, 



ii 
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were used. The scale is fixed by the phenomenological length ro [13] and the quark mass is 
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fixed by m^ro value. The parameters of our lattices are presented in the TABLE [H Note, 
that we used for T pc the value depending on m^ro which was obtained using interpolation 
of DIK results as function of rn^o 14|. For this reason values for T/T pc in TABLE |I] are 
slightly different from those in [g]. 

We evaluate a s following computational procedure used in 
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171 ] where pure gauge 



theories were investigated. Let us consider z direction as the 'time' direction of the spatial 
Wilson loops. The temperature dependent spatial static potential V S (R) is then defined as 

W(R, Z) 



aV s (R) = lim log 



(2) 



z^oo to W(R, Z + a)' 

where W(R, Z) are Wilson loops of size R/a x Z/a in lattice units. We assume an ansatz 
for V S (R) including linear and selfenergy terms: 



V s (R) = V + aR. 



(3) 



The ambiguous Coulomb term is not included and instead we make our fits at large distances 
R/a > 5 to avoid such term contribution. To increase the overlap with the respective ground 
state the APE smearing 18( of the links orthogonal to this direction was used. Furthermore, 
the hypercubic blocking (HYP) procedure 15!] was applied to links in ^-direction to decrease 
the selfenergy V and thus to decrease the statistical errors. We applied three steps of HYP. 
It occurred that the first step produced the maximal effect decreasing statistical errors of 
the potential V S (R) by factor about 5. Application of the second and third steps of HYP 
were much less effective but still useful, decreasing statistical errors by about 20%. The 
limit Z — > oo in (J2]) was approximated by taking the spatial Wilson loops with extension 
Z/a = 4 or 5 to evaluate V S (R). 



III. RESULTS 

In Fig. [T]we show the potential V S (R) with V s (ro) subtracted as function of R/tq for dif- 
ferent temperatures. One can see that potentials for all temperatures agrees with each other 
quite well, i.e. there is no strong temperature dependence in the temperature interval con- 
sidered. Respectively the spatial string tension is also virtually independent of temperature 
as can be seen from Fig. [2] In this figure we present our results for <y s {T) as a dimensionless 



ratio ^a s (T)/a, where a is the physical string tension at zero temperature [241 ] . We also 
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FIG. 1: Spatial static potential V S (R). 
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depicted in this figure DIK results for the finite temperature physical string tension cr(T) to 
emphasize drastic difference in the temperature dependence of cr(T) and cr s (T). 

Our results for a s (T) allow to conclude that for temperatures below T pc a s (T) is constant 
and equal to o. This does not come as a surprise since such equality was observed in the 
quenched theory long ago [lg, ll7j|. Unexpectedly the only previous result for T < T pc in 
un quenched lattice QCD 5| contradicts our result. In Ref. 5] yJa s (T) for T < T pc is roughly 
3T pc i.e. essentially higher then sqrta. The reason for such confusing result might be that 
in Ref. [s| different method to extract^the string tension a s was used. At T > T c our results 



are also different from those of Ref. [5j although agree qualitatively with results of Ref. [8|] 
Let us now compare our results with two theoretical predictions. In Ref. 21] the following 
expression for cr s (T) at T < 2T pc was derived within the field correlators approach 23]: 

a s (T) sinh(M/T) + M/T 



a 



cosh(M/T) + 1 



(4) 



where M is the inverse magnetic correlation length and is taken in Ref. 2lj equal to 1.5 GeV. 
Plotting square root of expression (jlj) in Fig. 2 we took M/T c = 7.5 thus taking T c = 
200 MeV. Note that in the considered temperature range the precise value of M/T c is not 
important. When M/T c » 1 equation (j4]) suggests very slow increase of the spatial string 



tension with temperature in nice agreement with our numerical results. Recently prediction 



for <7 S (T) was derived via AdS/QCD approach 



22|: 



°s{T) 



cxp 



1 



T > T r 



(5) 



and cr s (T) = a for T < T c . This expression also provide qualitatively good description of 
our data, as can be seen from Fig. 2. 




0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 

T/T pc 

FIG. 2: Square root of the ratio a s (T)/a (full circles) as function of temperature. The values of 
a are shown in Table I, see explanation in the text. Respective data for the finite temperature 
physical string tension a(T) (empty circles) computed by DIK 



curves are theoretical predictions from Ref. 



20| is shown for comparison. The 



2l| (thin curve) and Ref. 



221 ] (thick curve). 



We shall note that in our computations we do not keep the pion mass fixed, i.e. the 
computations are not made along the constant physics line. This may introduce sizable 
systematic effects if a s /a has strong dependence on the quark mass. In Ref. [5] a s was 
computed for two values of m n /mp, 0.65 and 0.85. The mild dependence on the quark 
mass was observed. These results need to be confirmed. In any case our conclusion about 
independence of a s on the temperature below transition is intact since m^ro varies by only 
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about 10% when temperature increases from T/T pc = 0.78 to T/T pc = 1.01, as can be seen 
from the Table [B 

Let us summarize our main results. We presented first numerical evidence in unquenched 
lattice QCD that for temperatures below T pc cr s (T) is constant and equal to a, the physical 
string tension at T = 0. This is in agreement with quenched lattice QCD results and 



2l| and eq© |22| for a s (T) at 



theoretical predictions. The theoretical predictions eq.([ 
temperatures T pc < T < lAT pc where all known lattice results clearly deviate from the 
expected high temperature behavior 



a s (T) = cg\T)T 



(6) 



are in good agreement with our results with better agreement for eq.( 
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